Comment on "Elasticity Model of a Supercoiled 
DNA Molecule" 

Bouchiat and Mezard 0] have proposed a "worm-like 
rod chain model" to describe recent experiments on DNA 
mechanics 0] . Their model can be solved using an anal- 
ogy with quantum mechanics, thanks to the local nature 
of the effective Hamiltonian. The letter concluded that 
this model has exceptional physical behavior: The con- 
tinuum limit is singular and linear response theory can 
not be used in the absence of a cut-off to calculate the 
experimental curves: For example the curves of extension 
versus torsion flatten in the continuum limit. 

The writhe, Wr = Xwl^i °f a cnam is given by the 
Calugareanu- White formula 0] 
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where r(s) = Jp t(s') ds' is the position in space and L is 
the length of the chain. Fuller showed that if we represent 
the tangent vector t by a point on the unit sphere the 
area, Xw^ enclosed by t(s) on the sphere Q, is in certain 
cases equal to Xw- The general relation between Fuller's 
and Calugareanu's formulas is expressed by: 


Xw = Xw mod. 47r - 


(2) 


Fuller's result, expressed as a local integral, was used in 
jjjj. The authors justify the use of Xw ^° interpret exper- 
iments on open chains with a topological modification of 
the sphere [pi and find that the probability distribution 

Pi.Xw) varies as 1 / (x w ) for large x w due to winding 
of trajectories about the south pole of the sphere. 

We have performed a simulations to directly compare 
the distribution of writhe implied by each formulation: 
We numerically generate a large equilibrated ensemble of 
scmiflexible chains of eight times the persistence length, 
made of 87V straight links with t(0) and t(L) parallel. 
This assures that the path of t(s) on the sphere is a 
closed loop. In order to use the Calugareanu- White for- 
mula, the chain must also be closed in real space. To do 
this we extend the chain at each end with a long straight 
segment in the directions — 1(0) and t(L) and then join 
the two extremities with an arc of a circle; the contribu- 
tion from the exterior is important in order to preserve 
the modulo equality eq. (§); it is this augmented version 
of eq. (1) which determins the rotation angle measured in 
experiments on open chains. We calculate Xw an0 - Xw 
for each chain from our ensemble and plot the integrated 
probability P(xw) = J*^p(x)dx- F° r the Cauchy dis- 
tribution, one thus obtains Xw ~ T^p f° r l ar ge values 
oi Xw- 

Results are displayed in fig. |l]. The curves of P as a 
function of Xw indeed depend on discretization. We have 
additionally verified that the curves of Xw sca led by the 
factor Vln ./V converge very accurately to a master curve 


as predicted in Q . On the other hand, Xw converges to a 
limiting curve as N grows. We find no evidence of diver- 
gence of Ym as ytrp nor evidence of scaling with ymiV. 
From fig. [l] we see that even for N small, the Fuller for- 
mulation strongly overestimates the statistical weight of 
large writhe fluctuations. 
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FIG. 1. Inset: The integrated writhe distribution func- 
tion P(xw) computed with two different formulations of the 
writhe. \w n °t computed for N = 2700 due to computer lim- 
itations. The main figure is a zoom of the top right corner of 
the inset. One can compare the continuous evolution of Xw 
as JV grows with the stability of Xw m the same limit. An 
asymptote in Xw ~ jTTp observed in this figure corresponds 
to the Cauchy distribution. The singularity of Xw ls clearly 
weaker than y^-p. Statistical fluctuations of the curves are 
smaller than their width. 

To conclude we believe that torsional fluctuations of a 
scmiflexible chains are not divergent when going to the 
continum limit. A DNA chain has a persistence length of 
order 50 times its radius; it is in a regime where eq. (|l|) 
has reached the continuum limit to a very good approx- 
imation. 
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